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Abstract 



' We generalize the algorithm for semi-linear parabolic PDEs in Henry-Labordere 

■ [11] to the non-Markovian case for a class of Backward SDEs (BSDEs). By simulating 

(-H ■ the branching process, the algorithm does not need any backward regression. To prove 

that the numerical algorithm converges to the solution of BSDEs, we use the notion 
of viscosity solution of path dependent PDEs introduced by Ekren, Keller, Touzi and 
Zhang [5] and extended in Ekren, Touzi and Zhang [6, 7]. 
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, path dependent PDEs. 

\o 

: 1 Introduction 

(N 

. Initially proposed by Pardoux and Peng [15], the theory of Backward Stochastic Dif- 

ferential Equation (BSDE) has been largely developed and has many applications in 
control theory, finance etc. In particular, BSDEs can be seen as providing a nonlinear 
^ I Feynman-Kac formula for semi-linear parabolic PDEs in the Markovian case, i.e. the 

I solution of a Markovian type BSDE can be given as the viscosity solution of a semi- 

linear PDE. We also remark that this connection has been extended recently to the 
non-Markovian case by Ekren et al. [5] with the notion of viscosity solution of path 
dependent PDEs (PPDEs). 

Numerical methods for BSDE have also been largely investigated since then. The 
classical numerical schemes for BSDEs are usually given as a backward iteration, where 
every step consists in estimating the conditional expectations, see e.g. Bouchard and 
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Touzi [2], Zhang [20]. Generally, we use the regression method to compute the con- 
ditional expectations, which is quite costly in practice and suffers from the curse of 
dimensionality. 

A new numerical algorithm for a class of semi-linear PDEs was proposed recently 
by Henry-Labordere [11], using an extension of branching process. It is a classical 
result that the branching diffusion process gives a probabilistic representation of the 
so-called KPP (Kolmogorov-Petrovskii-Piskunov) semi-linear PDE (see e.g. Watanabe 
[19], McKean [14]): 

dtu{t,x) + -D'^u{t,x) + /3i^^aku''{t,x) - u{t,x)j = (1.1) 

k=0 

with a terminal condition u(T,x) = ip{x), where is the Laplacian on W^, and 
(«fc)o<fc<no is a probability mass sequence, i.e. > and X^fcio^fc ~ '^^^ above 
semi-linear PDE (1.1) characterizes a branching Brownian motion, where every par- 
ticle in the system dies in an exponential time of parameter /3, and creates k i.i.d. 
descendants with probability a^. More precisely, let Nt denote the number of parti- 
cles alive at time T, and (Z^)j=i^... ^at^ denote the position of each particle, then up to 
integrability, the function 

v{t,x) := Et,, [Ufj^ V(4)] 

solves the above equation (1.1), where the subscript t,x means that the system is 
started at time t with one particle at position x. This connection is also extended for 
a larger class of nonlinear ity, typically n", a € [0, 2], with the super diffusion, for which 
we refer to Dynkin [4] and Etheridge [9]. Moreover, this representation allows to solve 
numerically the PDE (1.1) by simulating the corresponding branching process. 

When the coefficients are arbitrary in equation (1.1) and the Laplacian is 
replaced by an Ito operator Cq of the form 

£ou{t,x) := fi{t , x) ■ Du{t , x) + -aa'^{t,x) : D'^u{t,x), 

Henry-Labordere's [11] proposed to simulate a branching diffusion process with a prob- 
ability mass sequence {pk)k=o,--- ,M, and by counting the weight he obtained a so- 
called "marked" branching diffusion method. A sufficient condition for the convergence 
of the algorithm is provided in [11]. In particular, the algorithm does not need to use 
the regression method, which is one of the main advantages comparing to the BSDE 
method. 

For PDEs of the form (1.1), Rasulov, Raimov and Mascagni [21] introduced also a 
Monte-Carlo method using branching processes. Their method depends essentially on 
the representation of its solution by the fundamental solution of the heat equation. 

The main objective of this paper is to generalize the algorithm in [11] to the non- 
Markovian case for a class of decoupled Forward Backward SDEs (FBSDEs) whose 
generators can be represented as the sum of a power series, which can be formally 
approximated by polynomials. Although the polynomial generators are only locally 
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Lipschitz, the solutions may be uniformly bounded under appropriate conditions, and 
hence they can be considered as standard decoupled FBSDEs with Lipschitz genera- 
tors. 

We shall use the branching process to give a numerical algorithm, where the branch- 
ing process is constructed by countable number of independent Brownian motions and 
exponential random variables. To bring back the numerical solution to the BSDE 
context where there is only one Brownian motion with the Brownian filtration, we use 
the notion of viscosity solutions of path dependent PDEs introduced by Ekren, Keller, 
Touzi and Zhang [5] and next extended in Ekren, Touzi and Zhang [6, 7]. Namely, we 
shall prove that the numerical solution obtained by the branching diffusion is a viscos- 
ity solution to a corresponding semilinear PPDE, which admits also a representation 
by a decoupled FBSDE as illustrated in [5]. Then the numerical solution is the unique 
solution of the corresponding BSDE by the uniqueness of the solution to PPDEs. 

The rest of the paper is organized as follows. In Section 2, we consider a class of 
decoupled FBSDEs whose generators can be represented as a convergent power series. 
We then introduce a branching diffusion process, which gives a representation of the 
solution of the FBSDE with polynomial generator. In particular, such a representation 
induces a numerical algorithm for the class of FBSDEs using branching process. Then 
in Section 3, we complete the proof of the regularity property of the value function 
represented by branching processes. Next, we complete the proof of the main repre- 
sentation theorem in Section 4. For this purpose, we introduce in Section 4.1 a notion 
of viscosity solutions to a class of semilinear PPDE, where there is no non-linearity 
on the derivatives of the solution function, following Ekren, Touzi and Zhang [6, 7]. 
Finally, we illustrate the efficiency of our algorithm by some numerical examples in 
Section 5. 

2 A numerical algorithm for a class of BSDEs 

In this section, we shall consider a class of decoupled FBSDEs whose generators can be 
represented as a convergent power series, which can be approximated by polynomials. 
Then for FBSDEs with polynomial generators, we provide a representation of their 
solutions by branching diffusion processes. In particular, the representation induces 
a natural numerical algorithm for the class of FBSDEs by simulating the branching 
diffusion process. 

2.1 A class of decoupled FBSDEs 

Let Q° := {uj G C([0,r],M'^) : = O} be the canonical space of continuous paths with 
initial value 0, F*^ the canonical filtration and A'' := [0,T] x Q'^. For every {t,uj) G A*^, 
denote := supo<5<t Then the canonical process B{oj) = {Bt{oj) := ujt,0 < 

t < T} for all UJ & Q^, defines a Brownian motion under the Wiener measure Pq. 

Let : A° ^ R'^ and a : AO ^ S'^ be —progressively measurable processes. 
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Suppose further that for every < t < t' < T and uj,uj' € Q^, 

\fi{t,io)-n{t',uj')\ + \a{t,io)-a{t',io')\ < C{^/\t - t'\ + \\u;t^. - uj[,^,\\t) (2.1) 

for some constant C > 0, and aa^{t,uj) > cold for some constant cq > 0. We denote, 
for every (i,x) S A", by the solution of the following SDE under Pq: 

Xs = Xs, Vs <t and = xt + / ^{r,X.)dr + / a{r,X.)dBr, Vs > t.(2.2) 

Jt Jt 

For later uses, we provide an estimate on the SDE (2.2). 

Lemma 2.1. There is a constant C such that for every t & [0, T] anc? (ti, xi), (t2, X2) G 



< C(l + ||xi||2 + ||x2||?)(|tl-t2| + ||xi-X2||?). 



Proof. It follows by the same arguments as in Lemma 2 and Theorem 37 in Chapter 
V of Protter [16]. □ 

Suppose that ^/^ : il*' — )■ M is a non-zero, bounded Lipschitz continuous function, 
and F : (t, x, y) G A° X M — )• M is a function Lipschitz in y such that for every y, F{-,y) 
defined on A'' is F''— progressive. We consider the following BSDE: 

Yt = ^(°'°X) + ^ Fis, '^^''X.,Ys)ds- ZsdWs, n-a.s. (2.3) 

where the generator F has the following power series representation in y, locally in 
(t,x): 

00 

F(t,x,2/) := /3(j]afc(t,x)/ - 2/), (t,x)GAO, (2.4) 

k=0 

for some constant /3 > 0, and some sequence {ak)k>o of bounded scalar F"— progressive 
functions defined A''. We also assume that every is uniformly 1/2— Holder-continuous 
in t and Lipschitz-continuous in oj. 

Denoting by |.|o the L°°(A0) -norm |afc|o, we now formulate conditions on the power 
series 

4(s) := J^lofclos^ and £(s) := /3 [jV'lo %(s|^|o) - s] , ^ > 0, (2.5) 
fc>0 

so as to guarantee the existence and uniqueness of a solution of the backward SDE 
(2.3). 

Assumption 2.2. (i) The power series io has a radius of convergence < R < 00, i.e. 

£o{s) < 00 for \s\ < R and Iq{s) = 00 for \s\ > R. Moreover, the function I satisfies 
either one of the following conditions: 

(£2) or, £(1) > 0, and £~^{s)ds = T, for some constant s > 1. 
(ii) The terminal function satisfies < R. 
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Proposition 2.3. Let Assumption 2.2 hold true, then the BSDE (2.3) has a unique 
solution {Y,Z) such that supo<«7^ \Yt\ < R, ¥q— almost surely. 

Remark 2.4. When 1^ = 0, the function £ in (2.5) is not well defined. In order 
to provide a sufficient condition for the power series representation, we can consider 
the BSDE (2.3) with terminal condition Yt = £. Define the corresponding function 
Qeis) := f3[e~^£o{es) — s] . Suppose that for some e > 0, Assumption 2.2 holds true 
with the corresponding function qe, then by comparison result of BSDEs, the BSDE 
(2.3) admits still a unique solution (1", Z) such that Y is uniformly hounded. 

In preparation of the proof, let us consider first the ordinary differential equation 
(ODE) oi p{t) on interval [0,r]: 

p' = £{p), with initial condition p{Qi) = 1. (2.6) 

Lemma 2.5. Let |^|o < R, then the ODE (2.6) admits a unique hounded solution on 
the interval [0, T] if and only if Assumption 2.2 (i) holds true. Moreover, in this case, 
we have < p{t) < Vt G [0,T] for some Rq < R. 

Proof. First, since the function £ is Lipschitz on [0,L] for every L < ^y^, then it 
follows by Picard-Lindelof theorem (see e.g. Chapter 2 of Teschl [17]) that there is 
Tma.x > such that the ODE (2.6) admits a unique solution p on [0, Tmax) and that 
liuit^Tmi,^ \p{t)\ = j:^^ > 1- Further, we observe that ^(0) > 0, which implies that 
p{t) > on [0 

j^max)- Then it is enough to prove that Tmax ^ T- 
In the case that £{1) < 0, we have p{t) E [0,1] for every t £ [0,Tniax) and hence 
Tmax = oo > T. Otherwise, when £(1) > 0, we have p{t) > 1 for t G [0,Tmax) and it 
follows by (2.6) and the fact that p{R/\ip\o) = oo that J(^ds = Tmax- We hence 
deduce that T < Tmax by Assumption 2.2 (£2). □ 

Remark 2.6. The ODE (2.6) can he rewritten as 

p{t) = p(0)+ f £{p{s))ds. 
Jo 

Let ip{t) := p[t)\il:\Q, then under Assumption 2.2 we have 

ft °° 

e'^V(t) = ^(0) + / e^'p{j^\ak\w\s))ds. (2.7) 

fc=0 

In other words, the existence and uniqueness of solution to (2.6) is equivalent to that 
of (2.7). 

Remark 2.7. Suppose that = for every k > hq with some no G N, then 
clearly £{s) := p{Y.IU H\oH\o~^s^ - s) . Denote £,{s) := /3( Efcio 1(1 + ^)«fc|o|(l + 
e)'ip\Q~^ s'^ — s) . Let Assumption 2.2 hold true, then for e > small enough, we still 
have ^e(l) < 0, or ^^(1) > and J^^ J^^^ ~ ^ /^'^ some constant 1 < < oo. It 
follows that the ODE: p'{t) = £e{p) with initial condition p{Q) = 1 admits a unique 
solution on [0,T] under Assumption 2.2. 
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With the above existence and uniqueness result of ODE (2.6), we get the existence 
and uniqueness of the BSDE (2.3) in Proposition 2.3. 



Proof of Proposition 2.3. By Lemma 2.5, the solution p of ODE (2.6) is uniformly 
bounded by |^ with some constant C = Rq < R, where R is the convergence radius 
of the power series J2T=o I'^fclo^;^- Denote yc '■= —C V {y AC) for every y G M, 



^c(s,x,y) := F(s,x,yc) 

and 

oo oo 

lc{s,y^,y) ■= l3{^J2\ak\o\yc\'' - yc^, f_c{s,^,y) := -/^(^ |afc|o|ycl'' + |ycl 

fc=0 fc=0 

Then Fc, fc ^-nd / are all globally Lipschitz in y, and / < Fq < fc- Moreover, 
if we replace the generator F by /(--(resp. /„), and the terminal condition ^0 by IV'lo 
(resp. — iV'lo) in BSDE (2.3), the solution is given by Z := (resp. Z := 0) and 

Yt := p{T-tMo (resp.y, := -p(T-tMo). 

Therefore, by comparison principle, it follows that the solution (YcZq) of BSDE 
(2.3) with generator fc satisfies Y_ < Yc < Y, and hence lie*! < C. Further, since 
F{t,x,y) = Fc{t,x,y) for all \y\ < C, it follows that {Yc,Zc) is the required solution 
of BSDE (2.3). □ 



2.2 A branching diffusion process 

Let /3 > 0, no > and p = {pk)o<k<nQ be such that Y.k<noPk = 1 and pk > 0, k = 
0, • • • , no. We now construct a branching diffusion process as follows: a particle starts 
at time t, from position x, performs a diffusion process given by (2.2), dies after a 
mean /3 exponential time and produces k i.i.d. descendants with probability pk- Then 
the descendants go on to perform diffusion process defined by (2.2) driven by inde- 
pendent Brownian motions. Every descendant dies and reproduces i.i.d. descendants 
independently after independent exponential times, etc. In the following, we shall give 
a mathematical construction of this branching diffusion process in three steps. 



A birth-death process Suppose that we are given a probability space {Q, F, P), in 
which {T^'-')ij>o is a sequence of i.i.d. random variables drawn from the £{/3) exponen- 
tial distribution with mean /5 > 0, {In)n>i is a sequence of i.i.d random variables which 
are independent of {T^'-')i,j>o, of multi-nomial distribution A4(j>), i.e. P(Ii = k) = p^, 
Vfc = 0, 1, • • • , no. We shall construct a continuous time birth-death process associated 
with the coefficient (3 > and the probability density sequence (pfc)o<fc<no- 

The branching process starts with a particle at time 0, Nt denotes the number 
the particles in the system, every particle runs an independent exponential time and 
then branches into k i.i.d. particles with probability pk- We denote by T„ the n— th 
branching time of the whole system, at which one of the existing particles branches 
into In particles. Between T„ and Tn+i, every particle is indexed by (fei,-- - ,A;„) G 
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{1, • • • , rao}"', which means that its parent particle is indexed by (/ci, • • • , between 
Tn-i and T„. We also have a bijection c between N and U„>i 

2{i. -."or defined by 

n 

c((A;i,-- - := Y^U^o + lT- (2.8) 

Denote by ICt the collection of the indexes of all existing particles in the system at 
time t. Then the initial setting of the system is given by 

No = h To = o, ri = r0'0, ]Ct = {{!)}, yt£[o,n], 

and we have the induction relationship 

A^T,+i = NT^+Ii+,-l, r,+i = Ti+ min = + ^'^(^'+1), 

where Ki^i denote the index of the particle which branches at time Tj_|_i. Let 

/Cr,+i := {iK,+i,m) : 1 < m < U {(k,!) : k e JCtA {K^+l}} ■ 

In particular, if Ij+i = 0, then /CTi+i = {{k,l) : k G /Ct- \ {A'j+i}}. Clearly, at 
a branching time Tj, the particle branches into li particles which are indexed by 
{Ki, !),••• , {Ki, li), and all the other particles with index k are re-indexed by {k, 1). 
Let 

Nt := and ICt := /Ct^, for all t £ [Ti,Ti^i). 

Then {Nt)t>o is a continuous Markov process taking value in N. Since it is possible that 
a particle dies with k = descendants, the branching system is subject to extinction 
in finite time horizon, i.e. W[Nt = for some t > 0] > 0. Furthermore, {JCt)t>o is a 
random process taking value in □^6^2''"'^' "'""-^", and Nt = whenever ICt is empty. 

Example 2.8. We give an example of the branching birth- death process, with graphic 
illustration below, where hq = 2. The process starts with one particle indexed by (1), 
and branches at time Ti , • • • , T5 . The index of the branched particles are respectively 
(1), (1,1), (1,2,1), (1,1,2,1) and (1,1,1,1,1). At terminal time T, the number of 
particles alive are Nj- = 5 and 

ICt = {(1,1, 1,1, 1,1), (1,1, 1,1, 1,2), (1,1, 2, 1,1,1), (1,2, 1,1, 1,1), (1,2, 1,2, 1,1)}. 

• At time Ti, particle (1) branches into two particles (1,1) and (1,2). 

• At time T2, particle (1,1) branches into (1,1,1) and (1,1,2), particle (1,2) is 
reindexed by (1,2,1). 

• At time T3, particle (1, 2, 1) branches into (1, 2, 1, 1) and (1, 2, 1, 2), the other two 
particles are reindexed by (1,2,1). 

• At time T4, particle (1,1,2,1) branches into one particle (1,1,2,1,1), the other 
particles are reindexed. 
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• At time T^, particle (1, 1, 1, 1, 1) branches into (1, 1, 1, 1, 1, 1) and (1, 1, 1, 1, 1, 2), 
the other particles are reindexed by (1, 1, 2, 1, 1, 1), (1, 2, 1, 1, 1, 1), (1, 2, 1, 2, 1, 1). 



(1, 1, 1, 1, 1, 1) 




Lemma 2.9. For every probability density sequence {pk)o<k<no, we have lim„_j.oo Tn = 
00, a.s. In particular, the system is well defined from to 00. 

Proof. Without loss of generality, we can consider the case when pk = 0, VA; < uq and 
PnQ = 1. We first claim that Nf < 00 for all t > 0, it follows that sup{n : T„ < t} < 00 
for alH > and hence lim„^oo Tn = 00. Then to conclude, it is enough to prove that 
Nt < 00, \/t > 0, which means that the population of the particles never explodes. 
In fact, it is enough to use Example 2 of Kersting and Klebaner [12] to finish the 
proof. □ 



The branching Brownian motion Suppose that in the same probability space 
{Q,J-,¥), there is a sequence of independent d— dimensional standard Brownian mo- 
tions (W^,W'^,- ■ ■), which is also independent of the exponential random variables 
(r*'-^)jj>0 and multi-nomial random variables (/n)n>i- We can then construct a 
branching Brownian motion which starts at time t > 0. 

For the first particle in the system indexed hy k = (1), we associate it with a Brown- 
ian motion on [t, 00), defined by BlfJ = W^,yO<s< Ti. Let k = {h, ■ ■ ■ , kn) € /Cr„ 
be the index of a living particle at time T^, whose parent particle is indexed by 
(fci, • • • , kn-i), we associate it with a Brownian motion between [t, t + T^+i], defined 

by 

By the strong Markov property of the Brownian motion, it is clear that conditioned 
on (T*'-')jj>o and (/n)n>0; every process {Br^)r>t for k G )Ct is a Brownian motion. 
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In particular, given two particles = • • • , k^) and k'^ = {kf, • • • , k"^) such that 
kj = kj for all j = 1, • • • , i, the associated Brownian motion i?*''' and i?*''' share the 
same path before time t + Tj. 

The branching diffusion process To construct a branching diffusion process, 
we first remark that for every (t,x) G AO, the SDE (2.2) with initial condition = 
Xs , < s < t has a unique strong solution *'^X adapted to the natural Brownian filtra- 
tion, hence there is a progressively measurable function <1>*'-^ : [t, T] x C{[t, T], M*^) — t- M 
such that ^'""Xs = ¥'""{3, B.), Pq - a.s.. 

Then a branching diffusion process ^^^X^ is given by 

^'""Xj^^^ := + s,S.*'^'), VsGM+ and A; G /C^. (2.9) 

Moreover, for later uses, we extend ^'^X^^^ on the whole interval [0, T] by 

^'''XW := X, Vs < t and := ^>*'^(s, 5*'^^^), Vs > i. 

Remark 2.10. By the flow property of the SDE (2.2), we have that for every (t,x) G 
A^, r < s and k G JCs, 

¥'^{t + s,{B'/)u>t) = ¥+'-''-''''\t + s,iBl;%>t+r), r-a.s. (2.10) 

To conclude this subsection, we equip the above system with two filtrations. First, 
F = {Tt)t>o with 

J^t ■■= cr{{Tn,In,Kn)lT„<t + dlT„>t, n>l), 

where d denotes a cemetery point. Intuitively, F is the filtration generated by the birth- 
death process. In particular, T„ is a F— stopping time and Tt„ = a[{Tii,Ik, -K^A:)i</c<n) • 
Next, for every t > 0, let F* = {Tl_^g)s>o be the filtration on the probability space 
(ri, J^, P) generated by the branching diffusion process, which is defined by 

TUs ■■= V ^i^r, B^^\ B\'^^, 0<r<s, k€ K,). (2.11) 

2.3 Branching diffusion representation of backward SDE 

Using the branching diffusion process defined above, we can provide a representation 
of the solution to the decoupled FBSDE (2.3). 

Let (t,x) G A^, we consider the branching diffusion process {^'^X^)k£KT '^^ [^i^] 
defined in (2.9), where the probability sequence p = {pk)o<k<no satisfies that p^ > 
whenever |afc|o 7^ 0. Denote 

Wt,x :=n;Lr( "'"^^^^"' ''-X^-) y ^^^^^ MT-*:=sup{n:t + r„<T} 

is the number of branchings occurred in the particles system before time T. Our main 
representation formula is the following function: 

v{t,^) := ¥f [ft,.] with ^t,. := Wt,.UkeKT-tH''^X''), (2-12) 
where the integr ability of ^t,yc is verified in the following result. 
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Proposition 2.11. Suppose that Assumption 2.2 holds true. Then for every (i,x) G 
A*^, the random variable ^^.x given in (2.12) is integrable and the value function v is 
uniformly bounded. Moreover, for every M > 0, there is a constant C such that 

\v{t,uj) - v{t' ,u:')\ < C(V|i-t1 + ||wtA. -^I/a-IIt), 

whenever \ {t,uj)\ < M and \ {t',uj')\ < M. 

The proof of Proposition 2.11 will be completed later in Section 3. 
Our main result of the paper is the following representation theorem. Let ^'^X be 
the unique strong solution to the SDE (2.2) in the probability space {Q,T,F), denote 

y^o ^(t^x). (2.13) 

We also consider the BSDE (2.3) with generator 

no 

F,,,{t,^,y) := /3(J^afe(t,x)/ - yj. (2.14) 
fc=0 

We define £no by 

no 

^no(s) := (3[^\ak\om'^-h'-s), Vs > 0. 

fc=0 

It is clear that when Assumption 2.2 holds true for i, then satisfies also Assumption 
2.2. It follows from Proposition 2.3 that the BSDE (2.3) with generator Fn^ has a 
unique solution, denoted by (Y, Z), such that Y is uniformly bounded. 

Theorem 2.12. Suppose that Assumption 2.2 holds true, and {Y, Z) is the unique 
solution of BSDE (2.3) with generator Fn^ (defined by (2.14) ) such that Y is uniformly 
bounded by Rq, the constant introduced in Lemma 2.5. Then Y^ = Y, FQ — a.s. 

The proof of this result will be provided in Section 4 using the notion of viscosity 
solutions to a path dependent PDE. 

Remark 2.13. The results in Proposition 2.11 and Theorem 2.12 hold true for any 
probability sequence p = {pk)o<k<no satisfying that pk > whenever \ak\o 7^ 0. This 
implies that the integrability and expectation of'^t,^ is independent of the choice of p. 
However, the variance of ^t,x does dependent on p, where an upper bound is given by 

U^:-^(\^I^)u,,^,_,ml (2.15) 

^ Pin ' 

Using Proposition 2.11 and Theorem 2.12 with a simple manipulation on the coeffi- 
cients, and in spirit of Lemma 2.5 as well as Proposition 2.3, the upper bound (2.15) 
is finite if and only if £""{1) < or £^{1) > and {■s))~^ ds = T for some s > 1, 
where 

no I 1 2 

ris) :=/5(^^HrV-.). 

fc=0 ^'^ 
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2.4 Numerical algorithm by branching process 

The representation result in Theorem 2.12 induces immediately a numerical algorithm 
for BSDE (2.3) by simulating the branching diffusion process. For numerical implemen- 
tation, the branching times can be simulated exactly since they follow the exponential 
law, and the diffusion process can be simulated by a Euler scheme. 

Let A = {to, ■ ■ ■ ,tn) be a discretization of the interval [0,T], i.e. = to < • • • < 
tn = T. Denote |A| := maxi<fc<„(tfc — tfc„i). To give the Euler scheme, we introduce 
the frozen coefficients ji^ and by 

^^(t,x) := i^L{tk,^^) and cT^(t,x) := a{tk,±^), Vt G [t^, t^+i), 

where x^ denotes the linear interpolation of (x^g, • • • , xt„) on the interval [0, T\. Then 
clearly the process given by the SDE 

X^ = [ ^i^{s,X^)ds + [ a^{s,X^)dBs, Po- a.s. (2.16) 
Jo Jo 

can be simulated, which is also the Euler scheme of the SDE (2.2). By standard 

arguments using Gronwall's Lemma (see e.g. Kloeden Platen [13] or Graham and 

Talay [10] in the Markov case), we have the following error analysis result: Let X be 

the solution process of (2.2) with initial condition (t,x) = (0,0), X^ be the solution 

of (2.16) and X^ denotes the linear interpolation of {Xf^, ■ ■ ■ ,Xf^) on [0, T]. 

Lemma 2.14. There is a constant C independent of the discretization A such that 



E 



sup l\Xt-Xf^\' + \X,'^-Xf^f 



< C lAI. 



-0<t<T 

Moreover, for the BSDE (2.3) with a general generator function F : [0, T] x il*^ x M — )• 
M which admits a representation (2.4), we can approximate it by some polynomial Fn^ 
of the form (2.14). Let F^^{t, x, y) := /?( ZZo «^(*' " v) ' ^^ere 

a^(t,x) := afc(ti,x^) for every A; = 0, •••,no and t£[ti,ti+i). 

Further, under Assumption 2.2, by simulating the branching diffusion process (X^''^)^^^^ , 
the numerical solution 



„ A (rp ^AjKn \ 



Yrf" := E 

L PL 
is the solution of the following BSDE 

Yo = ^^^{X^) + r F^^{t,X^,Yt)dt - r ZtdBt, Po - a.s. 
Jo Jo 

Finally, we provide an error estimation of the numerical solution: 

Proposition 2.15. Under Assumption 2.2, there is a constant C independent of uq 
and A such that 

\Yo^-Yo\ < C(|F-F„Jioo(AOx[_^„,fi„]) + VA). 

Proof. This estimate follows from a direct application of the stability result of back- 
ward SDEs together with the error estimation in Lemma 2.14, see Proposition 2.1 and 
their subsequent remark in El Karoui, Peng and Quenez [8]. □ 
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3 Holder and Lipschitz regularity of v 



This section is devoted to the proof of Proposition 2.11. We first derive some estimates 
of the birth-death process defined in Section 2.2, then together with the tower property, 
we can complete the proof of Proposition 2.11. 



3.1 Some estimates of the birth-death process 



We recall that F = {J-t)o<t<T is the filtration generated by the birth-death process 
defined in the end of Section 2.2, and that the number of branchings occurred in the 
system before time t is denote by Mt := sup |n : Tn < t}. We also introduce: 



r]{t) := ¥F 
Lemma 3.1. For every < s < t, 



n 



n=l 



Pin 



n 



Mt \"-lJo 
Pin 



n 



M, |a/n|0 



n=l 



Pin 



and 



(nfir 



Pin 



|QJilo 
Ph 



{r]it-Ti)YnT,<t. 



Proof. Let Z he a random variable and A £ J^, then C^{Z) denotes the law of Z and 
C^[Z\A) denotes the distribution of Z conditioned on A under the probability P. We 
notice that for every i,j > 1 and s > 0, 



Let < s < i, the law of number of branches between s and t (which equals to Mt — Ms) 
is completely determined by A^^, {T^'-^)i>Ms,j>Q and {Ii)i>Ms+i- It follows that 

C^{Mt - Ms, {lMs+i)i>i \Ns = l, ICs, Ms = j) 
= C^{Mt - Ms, {lMs+i)i>i \Ns = l, ke JCs, Ms = j, T^^--^^^^) > s) 
= C^{Mt-s, {Ii)i>i), 

and hence 

C^{Mt - Ms, {lMs+i)i>i \Ns = l) = C^{Mt-s, {Ii)i>i). 
Since Nt = Ns + X]n=A/^+i(-^n — l)i we deduce that 



E^ 



TjMt 



Pin 



Ns 



r](t - s). 



Moreover, since every particle branches independently to each other, we deduce that 



E^ 



Y[Mt 



Pin 



1^' 



Ns 



ivit-s)y, 
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which imphes that 



n 



n=Ms+l' 



Pin 



And hence 



n 



Mt |QJn|0 
" = ' Pin 



E" 



n 



Ms |Q/n|0 
" = ' PI^ 



TTMt 



1^' 



For the second equahty, we remark that (/j)i>2 and (T^'^)i>2, j>o are all indepen- 
dent of (Ti, Ji). Then consider a family of conditional probabilities (lPs,j)seR+,ig{o,--- ,no} 
of P w.r.t. the a— field generated by (ri,/i). Under a probability P^^j, the law of 
{Mt,Nt)ls<t depends only on {Ij)j>2 and (r-''')j>2, />o- Considering in particular 
i = 1, we have 



And hence 



E^ 



ttA/j 



Pin 



S<t 



s<t- 



Moreover, by the independence of the evolution of i particles under Pg^j, we get 



E^ 



n 



Mt 

n=M a+l 



|a/„|o 



Pin 



{7]{t - s))h,<t. 



which concludes the proof. 



□ 



Lemma 3.2. Suppose that for some t > 0, r][t) < oo. Then there is 6 > such that 
r]{s) < oo, for every s £ [t,t + 6]. 

Proof. First, it follows from Lemma 3.1 that for every t, S > 0, 



ri{t + 5) 



E^ 



n 



Ms |QJn|0 



n=l 



Pin 



Ns 



Let us consider another pure birth process {Nt,ICt) on a probability space {(l,JF,F) 
with the same constant characteristics [3 and {pk)o<k<no such that p„p = 1. We suppose 
without loss of generality that uq > 2 and denote C := maxo<fc<no + ^(0- Then 
clearly it is enough to prove that E'^[C^*] < oo for some (5 > to conclude the proof. 
The distribution of Ns can be computed explicitly (see e.g. Athreya and Ney [1]) and 
verifies that for some constant C > 0, 

F[Ns = n\ < Ck"s with Ks := (i _ e-^/3("<'-i))^/("»-'). 

Then for 6 > small enough, ks is small enough such that E'^[C^'^] < oo. □ 
The birth-death system is closely related to the ODE (2.6). Let us define 



E^ 



(ivMi)(iviVj)(n^ii^ 



(3.1) 
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Proposition 3.3. Suppose that Assumption 2.2 holds true. Then 



sup r]{t) < oo and sup Dt < oo. 

0<t<T 0<t<T 



(3.2) 



Proof. We first observe that r/(0) = |V'|o by its definition, and it follows from Lemma 
3.1 that 

no 

77(t) = E[|^|olTi>t] + E[Y,\ak\o{v{t-Ti))''lT,<t 

k=0 
ft , "0 



k=0 

7._n 



fe=0 

Suppose that Tq := inf{s : r]{s) = oo} < T, then it follows from Lemma 2.5 and 
Remark 2.6 that 7?(t) = /9(t)|^|o < oo, where p is the unique solution to the ODE 
(2.6). This contradicts Lemma 3.2. We then have Tq > T and 'r]{T) < oo. Since r] is 
increasing, this provides the first claim in (3.2). 
We next denote 



%(t) := 



n 



Mt l(l+e)a/„|o 



n=l" 



Pin 



1(1 



Nt 



In spirit of Remark 2.7, we know that for e > small enough, r]f,[T) < oo. It follows 
that 



sup Dt := sup E 

0<t<T 0<t<T 



(1 vMt)(i viVt) n„, T, 



,<t- 



< OO, 



since there is some constant > such that n < Ce{^ + e)" for every n > 0. And we 
hence conclude the proof. □ 



3.2 Proof of Proposition 2.11 



In preparation of the proof, we first provide a tower property of the branching diffusion 
process. Let (t, x) G A'^ and r : 0'^ — )■ be a F''— stopping time such that t > t, then 
f := r( *'-^X.^^^) is a F*— stopping time in the probability space (fi, P). 

Lemma 3.4. Suppose that Assumption 2.2 holds true, let (t,x) G A*^, < s < T — t 
and f he given above. Then we have 



and 



v{t,x) 



n 



n=\ 



n 



,v{t + s, 



t,xj^k 



(3.3) 



V{f, ''^X^'^)lt+T,>, 



+ 



(t + Ti, *'-X«)lt+T,<f 



Ph 



(3.4) 
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Proof. First, following the arguments of Lemma 3.1, we know 

C^{Mt - Ms, (/M.+i)j>i, {W^'^+')i>i \Ns = h Ms = i) 
= C^iMt-s, {Ij)j>i, iW%>i). 
Together with the flow property of SDE in (2.10), it follows that 



=Ms+l 



Pin 



Ns 



1, Ms = i, k G ICs, C'^X^)t<r<t+ 



v{t + s, 



Then by the independence of evolution of every particle in /C^, (3.3) holds true. 

For the second equality, we consider a regular conditional probability distribution 
(r.c.p.d.) (Ptj)tjGno of¥ w.r.t. a{B^^^,) (see also Stroock Varadhan [18] for the notion 
of r.c.p.d.). Then using the fact that under P^, the law of Ti conditioned on Ti > f — t 
is the same as that of Ti , we get that 

E''"[*t,xlt+Ti>fH I i + Ti >f(L.)] = i;(f(a;),X.*'"'(')). 

Further, using similar arguments as in Lemma 3.1, by considering the distribution of 
^t,xlt+Ti<-f conditioned on -Fy^, we get 



IE^[^i,xlt+T,<f] = 
which concludes the proof. 



Ph 



□ 



Proof of Proposition 2.11. (i) First, it follows immediately from Proposition 3.3 
that ^'j^x is integrable and |f(t,x)| < p{T — t)|^|o < Ro- 

(ii) Let t £ [0, T] and xi, X2 € 0,^. By Lemma 2.1, for every s G [t, T] and k £ /Cg, we 
have 



sup 

t<r<s 



< C(l + ||xi||t + ||x2||j) ||xi-X2||t, 



for some constant C independent of xi, X2. Then using the fact that {ak)o<k<no and 
ijj are all Lipschitz in x, 

\v{t,xi) - v{t,^2)\ < E^[\^t,.,-%,^,\] < C AE^[|| 

< C(l + ||xi||t + ||x2||t) ||xi-X2||t, 

where Dt is defined in (3.1). 

(iii) Let < s < t < T, then it follows from Lemma 3.4 that 



\v{s,x) -v{t,XsA.)\ < 



E^ 



n 



Mt. 

71=1 



Pin 



< C{ sup Dr)E^[ sup Ixs-^'^'Xl^l 

s<r<t rels,t\ 



+ 



E^ 



^m-.ajjt^ (.(t,x.,.)) 



Nt-s 



Pin 



v{t,x 



< c{i + Ms)Vt^ + !</>(*) 
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where (p is the unique solution of the ODE 



no 

afc(t, x)i?i>'^(r) — 0(r) j with terminal condition (l){t) = v{t,yis^.). 

Moreover, by comparison principle of ODE, \(j){r)\ < p{r), Vr G [s,t]. Then \(p{t) — 
^ C{t — s) for some constant C independent of {s,t,x), which implies that v is 
locally (1/2) -Holder in t. □ 

Remark 3.5. When {ak)o<k<no O'^d if) are all constants, the value function v{t,x) is 
independent ofx andt i-)- w(r— x)|^|q ^ is a solution of the ODE (2.6). Therefore, in 
spirit of Lemma 2.5, Assumption 2.2 is also a necessary condition for the integrability 



4 The branching diffusion representation result 

This section is devoted to the proof of Theorem 2.12. 

We first consider a class of semi-linear parabolic path-dependent PDEs (PPDEs) 
and introduce a notion of viscosity solution, following Ekren, Keller, Touzi and Zhang 
[5] and Ekren, Touzi and Zhang [6, 7]. Our objective is to show that the value function 
V, defined by our branching diffusion representation, and the Y— component of the 
BSDE are viscosity solutions of the same path-dependent PDE. Then, our main result 
follows from a uniqueness argument. 



4.1 Viscosity solutions of PPDEs and FBSDEs 

We consider a PPDE which is linear in the first and second order derivatives of the 
solution function. This is a simpler context than that of [5, 6, 7]. As a consequence, 
following Remark 3.9 in [6], we use a simpler definition of viscosity solutions. We 
shall also provide an (easy) adaptation of the arguments in [6] which relaxes their 
boundedness conditions, thus allowing the terminal condition and the generator to 
have linear growth. 



4.1.1 Differentiability on the canonical space 

For ah t G [0,r], we denote by 0* := {to G C{[t,T],R'^) : ut = 0} the shifted canon- 
ical space, the shifted canonical process on 0*, F* the shifted canonical filtration 
generated by B^, Pq the Wiener measure on il* and A* := [t,T] x 0*. 

For s < t, u £ ^l^ and u' £ Q*', define the concatenation path uj ®t'^' ^ by 

(CJ w')(r) := UJrls<r<t + {uJt + ^^'r)U<r<T , VrG[s,r]. 

Let ^ G J-q^ and V he a F*^— progressive process, then for every {t,bj) G A*^, we define 
e*'^ G F't and {yt"')t<s<T by 
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Following Ekren et al. [6, 7], we define some classes of processes in A*, t > 0. Let 
C°(A*) be the collection of all F* —progressive processes which are continuous under 
the norm doo, where 

doo{{s,Uj),{s',u')) := \s-s'\ + sup jn^sAr - W^'Arl' V(s, w), (s', w') ^ A*. 

t<r<T 

Denote by C^(A*)(resp. C/C(A*)) the collection of functions in C'^(A*) which are uni- 
formly bounded (resp. uniformly continuous), and f/Cb(A*) := C/C(A*) n C°(A*). 
Next, denote by the solution of the SDE on (J]*, J^,P^): 

Xs = X,, Vs <t and Xs = xt+ [ ii{r,X.)dr + [ a{r,X.)dBl, Vs > t. (4.1) 

Jt Jt 

Clearly, under has the same law as that of introduced in (2.2) under Pq. 

We denote the induced measure on the shifted space 0* by: 

Pi,x:=Poo(xO'*'''-xt)"' and Px := Po,o. (4.2) 

Remark 4.1. Let {t, x) G A'', t > t be a ¥^ — stopping time on il*, ^ G ^'^^ (Pi^)LJGf7 
6e a regular conditional probability distribution (r.c.p.d., see Stroock-Varadhan [18]) of 
Pi,x w.r.t. then clearly, = e'^-(-)>" [^^H'"^] forWt^^-a.s. a; G J]. 

For every s G [0, T) and u : — > M, we introduce the Dupire [3] right time- 
derivative of u defined by the following limit, if exists, 

» r+ \ r u{t + h,uj.At)-u{t,uj) a ^ 

Otuit, uj ):= hm ; , t<i, and Ojuf i , w) := hm Otu[t,uj). 

hiO h t<T,t^T 

Definition 4.2. Let u be a process C°(A*). We say u G C^'^(A*) if dtu G C°(A*) and 
there exist d^u G C°(A*,M'^), d^^u G C°(A*,§'^) suc/i i/iai, for any {s,uj) G A*.- 

dul^^ = {dtu)P'^dr + {d^u)P'' ■dB'^ + ^{d^^u)P'' ■.d{B')r, Pt,x-a.s. (4.3) 

If, in addition, u G C^(A*), we i/ien say u G C^'^(A*). 

It is clear, for s < t, w G 0° and n G C^'^(A''), we have u*''^ G C^'^(A*). 

Finally, for all t G [0, T], we denote by T* the collection of all F*— stopping times r 
such that {uj : t{uj) > s} is an open set in (fi*, || • H^) for all s G [t, T], and by 7^ the 
collection of stopping times t £ T*" such that t > t. The set A*(t) := {(t,(^) G A* : 
t < t{uj)} is the corresponding localized canonical space, and we define similarly the 
spaces CO(A*(t)), Ci'2(A*(r)), etc. 

4.1.2 A path-dependent PDE 

In this section, we do not need the restriction of the generator to have a power series 
representation in y as in (2.4). We then consider a slightly more general generator 
i*" : A'^ X R — )• M such that {t,oj) i — ?■ F{t,(jj,y) is F*^— progressive for every y G M. 
Consider the second order path-dependent differential operator: 

Cif := ^ • d^if + : dl^ip. (4.4) 
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Given a J^j'— measurable r.v. ^ : — t- M, we consider the path-dependent PDE: 

{-Cu- F{.,u)}{t,u;) = 0, t<T,uj£n^, (4.5) 
u{T,uj) = wGJ]°. (4.6) 

Assumption 4.3. There is a constant C such that supj<2^ |F(t,0,0)| < C, and 

\F{t,u:,y)-F{t,J,y')\ + \i{uj)-i{J)\ < C{\y - y'\ + - oj'\\t) , 

for every t e [0,T], {uj,y), {J ,y') G 0° x M. 

We next introduce two classes of test functions for every F'^— adapted process w. 

Au{t,uj) := {v? G C7i'2(Ai) : 3 H G r^, (99 - n*'-)t(0) = mm WF^'- [{if - n*''^),AH] }, 
Au{t,uj) := |(^ G Ci'2(A*) : 3 H G r_f, - u*''^)t(0) = niaxE^'-" [(99 - n*'^)^AH] }• 

The next definition requires the following notation for the path-dependent second order 
differential operator on the shifted canonical space: 

(£*'-^)(s,a;') := {^'''^ ■ d,^){s,u:') + ^{{aa^Y'^ : dlip){s,co'), {s,u;')€A'. 

Definition 4.4. Let u : A'^ — > M &e a locally bounded ¥^ — progressive process. 

(i) We say that u is a viscosity subsolution (resp. supers olution) of PPDE (4.5) if, for 
any {t,tj) G [0, T) x fl^ and any if G Au{t,uj) (resp. ip G Au{t,ijj)), it holds that 

{ -atV9-£*'^(/j-F*'^(-,n*'^)}(t,0) < (resp. >) 0. 

(ii) We say that u is a viscosity solution of PPDE (4.5) if it is both a viscosity subso- 
lution and a viscosity supersolution. 

Remark 4.5. (i) In Definition 4-4j can only use, without loss of generality, the 
test functions (p £ A (resp. A) such that (ip — u*'^)t{0) = 0. 

(ii) Similar to Remark 3.9 of Ekren, Keller, Touzi and Zhang [5], we can easily verify 
that under Assumption 4-3, for every X & W, u is a viscosity solution to (4.5) if and 
only if u{t,uj) := e'^'^'^~'^^u{t,uj) is a viscosity solution of 

-dtu- Cu- Fx{.,u) = 0, where Fx{t,uj,y) := -Xy + e^^F{t,uj,e~^^y) . 

(iii) Similar to Remark 2.11 of [5], we point out that in the Markovian setting, where 
the PPDE (4.5) decreases to a classical PDE, a viscosity solution in sense of Definition 
4.4 is consistent to the viscosity solution in standard sense, by the uniqueness result 
proved below. 

4.1.3 The existence and uniqueness of solutions to PPDE 

This section follows closely the arguments of [5, 6, 7]. However, their results do not 
apply to our context, because of the possible unboundedness of fi and a. Moreover, 
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the PPDE in our context is linear in the gradient and the Hessian components, which 
significantly simplifies the approach, see Remark 3.9 of [6]. 

The above viscosity solution to PPDE (4.5) is closely related to the following de- 
coupled FBSDE. For every (t,x) G A°, let X^'^'"" be the solution of (4.1), (yO'*'^, Z^'*'^) 
be the solution of the BSDE on (0*, J4,P*), 

Ys = e(X°'*''')+ / F(r,X.°'*'^,y^)dr- / Z^. ■ dB^ (4.7) 

J s J s 

By the Blumenthal 0-1 law, Y^'^'^ is a constant and we then define 

n(t,x) := y,°'*'^. (4.8) 

For later use, we observe that, since the diffusion matrix a is nondegenerate, the above 
BSDE (4.7) is equivalent to the following BSDE on (0*, J"*,, Pj,x): 

= e^^iB') + f F'^^{r, B'X)dr - f Z^ ■ {dBl - ^*'"(r, B')dr) , 

J s J s 

whose solutions are denoted by (Y^'^'^, Z^'^'^), i.e. yI^'^'^ = Y^'^'^ = n(t, x) for every 
(t,x) G A''. Moreover, by equation (4.6) of [6], we have the dynamic programming 
principle 

y^o,i,x ^ ui'^(T,B*) + r F^'''{r,B\Y^°'^''')dr 

J s 

- J^Z^'''^ ■ {dBl - fi''^{r, B')dr) , Ft,x - a.s., (4.9) 

for all (t,x) G A° and r G 7^. 

Theorem 4.6. Let Assumption 4-3 hold true. 

(i) There is a constant C > such that \/{t,uj), {t',uj') G A*^, 

\u{t,u) -u{t',uj')\ < C{\\uj\\t + \\uj'\\t'){^/\t-t'\ + \\uJtA- -^'fA-Wr)- 

(ii) u is a viscosity solution to PPDE (4.5). 

Proof, (i) is proved in Proposition 4.5 of Ekren et al. [6], since our BSDE (4.7) is 
a particular case to their equation (4.4). It is in fact an immediate consequence of 
Proposition 2.1 in El Karoui, Peng and Quenez [8] together with the estimation in our 
Lemma 2.1. 

(ii) We adapt the arguments in the proof of Proposition 4.5 [6] to our context. We 
only show that u is a viscosity subsolution. Assume u is not a viscosity subsolution, 
then there exist (t,w) G A° and ip G Au{t,uj) such that c := ^C*'"^ (/?(*, 0) > 0. Without 
loss of generality, we may also assume that ip{t,0) = u{t,u}). Denote, for s G [t,T], 

Y::=^{s,B'), Z',:=dMs,B'), 5Y, := Y^ - Y,'^'''^ , dZ^ := Z', - Z^/'^ . 

Applying Ito's formula together with Definition 4.2, we have for some F*— progressive 
process a bounded by some constant Lq, 

d{6Ys) = - [(£*'^vj)(s, S*) + as5Y} ds + • (dS* - ;U*''^(s, B^:)ds) , P^,^ - a.s. 
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Observing that 5Yt = 0, we define 

H := TAmi{s>t:C^''^ip{s,B^)-Lo\6Ys\<c/2}. 
Then by the continuity of 5Y as well as C^'^ip, H G 7I[; and 

5*) + as6Ys > c/2, for all s £ [t,H]. 
Now for any t £ such that r < H, we have 



= 6Y, 



6Yr + 



> ip{T,B^)-u^^'^{T,B^) + c{T-t)/2- 6Zs- {dBl- ^t^''^{s,B^)ds), 

'Pt,uj — a.s., where we also used the dynamic programming principle in (4.9). Then 
lEf*'" [{^ - 'u*''^)(t, B^)] < 0, contradicting the fact that ip G Au{t, to). □ 

For the proof of our main result, we also need the following comparison result. 

Theorem 4.7. Let Assumption 4-3 hold true, u^,u^ be two ¥— progressive cadlag 
processes on with corresponding jumps ssume that (resp. u^) 

is a viscosity subsolution (resp. supers olution) of PPDE (4.5), and u^{T, ■) < ^(•) < 
u^{T, •). Then < on A°. 

In preparation of the comparison principle, we first introduce two extended spaces 
'^(A*) and C< '^(A*) of C^'^{A°) and derive a partial comparison principle as in [6, 7]. 

Definition 4.8. Let (t, w) G A°, n : A* M be ¥*- adapted. 

(i) We say u G Cj'^(A*) if there exist an increasing sequence of ¥^ — stopping times in 
: t = tq < Ti < ■ ■ ■ < T such that, 

a) Ti < Tj+i whenever Ti < T, and for all Co G fi*, the set {i : Tj(a)) < T} is finite; 

b) For each i > and Co e n\ T-^f^'^ G T^^'-^^ and u^d^^'^ G Cl'"^ {A^^'-^^T-^^f'^''^)) ; 

c) u has non-negative jumps (Au >0), and 



i>0 



dtu + a^'^l + y^^d^ul ){s,B')ds 



< oo. 



(4.10) 



(fi) We say u G C;;:(A*) if -u G C:,';'(A*). 

Lemma 4.9. Suppose that Assumption 4-3 holds true. Let be a viscosity subsolution 
and u? be a viscosity supersolution of PPDE (4.5) such that u^{T,-) < it^(T, •). If 
e Cli{A°) or n2 G cJ;o(A°), then < on A". 

Proof. We follow the lines of Proposition 4.1 of Ekren, Touzi and Zhang [7]. Suppose 
that G CiQ'g(A'^). First, let us show that, for every i > and oj G fl*^, 



( 



(4.11) 
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Without loss of generality, it is enough to consider the case i = 0, where Pro(Lj),a; = 
for all G Assume to the contrary that 

2Tc := iul-uir - IE"^[«-<)1 > 0, 

we set 

Xt := (uj - ul)+ + ct, Yt := supper' Ei[^r An], t* := mi{t > : Xt = Yt} < n, 

where Et[C]{uj) := E^^'^iC^'"^] = E^x[C\Tt]{io). In particular, Eo[-] = E^^[-]. Then 
{yt)t>o is a supermartingale, {Yt/\r*)t>o is a martingale and r* is an optimal stopping 
time for the problem sup^g^-o Eo[Xt-]. It follows that 

Eo[X,.] = Eo[y..] = Yo > Xo = 2rc + Et[(ni^ > rc + Eo[X,J. 

Then there exists uj* € fl^ such that t* := t*{uj*) < ri. And therefore 

{u^ - uY{t*,u;*) + ct* = Xt^ioj*) = Yt,{uj*) > Et* [X^J > ct*, 

which implies that < {u^ - u'^)^{t* ,uj*). Set ip{t,uj) := (n^)**''^* (t, tj) + c(t*). Then 
if G C^'2(A**(ti)) since £ C^'2(A(ri)). Moreover, let 

H := inf {t>t* : u] - < O} A n G Tf . 

Then for every t £ , 

(^-(^z^f '-*)(r,0) = Xt,{u;*) > Kt'[YrAn]{u*) 

> E,.[Xr,,,]iu*) = E^**,-[(y,-(nY'"*).J, 

which implies that ip G Av?{t*,uj*). It follows that 

< -Cip{t*,uj*) < - c - Cu^{t*,co*), 

which contradicts the fact that is a subsolution and we hence prove (4.11). Further, 
since (lPri(a;),Lj)wer2 induces a r.c.p.d. of w.r.t. Jv- (see Remark 4.1), it follows by 
(4.11) that for every i > 0, 

By sending i — )• oo, we get that (n^ — m^)o < Eo[(u"'^ — u^)^] = 0, which completes the 
proof of Mq < Uq. □ 

Proof of Theorem 4.7. We follows the lines of the proof of Theorem 7.4 of Ekren, 
Touzi and Zhang [6], where a comparison principle for PPDE (4.5) was proved in case 
a = Id- In spirit of Remark 4.5, we suppose without loss of generality that F decreases 
in y. 

For every e > 0, we denote 
Oe := {x GR'^ : \x\ < e}, Oe ■= {xGM'^ : \x\ < e}, dOe := {x G : \x\ = e}; 
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Ol := [t,T) X Oe, Ol := %T] x 0„ dOl := {[t,T] x 90^) U ({T} x O,). 

Let to = 0, xo = 0, (ti)i>i an increasing sequence in (0,T] with tj = T when i is large 
enough, and (xi)i>i a sequence in M''. Set tt := (ti,Xi)i>o and Tin '■= {ti,Xi)Q<i<n- 
Given vr^ and (t, x) E Of^^, define 

Ho'"^'^ := inf{s > t : |B* + x| = e} A T, H^l^^'f := inf{s > H*'"^'^ : |S* - = e} A T. 

For t E T], let B^''^"'^'^{uj) denote the linear interpolation of {ti, X^j=o Xj)o<i<n and 
(h-'"''M,E"=o^j + ^ + ^H'--.^H)^>g- Define 

e^(7r„;(t,x)) := 3^,^'-"'*■^ 
where, omitting the superscripts ^'"^^'^'^^ y is defined 

ys = i{B)+ I F(r, ^.^^*,.,a[jjM,._H*-i^),x)dr- / Zr{dBr - f^ir, B)dr) , 

with H^^'^ := t. Then clearly, for every n and -7r„, the deterministic function 6"^ : = 
On^TTn, •) is the viscosity solution of the standard PDE on Of^: 

-dtOf, - fi{s,u;^-)D9f^ - ^aa^{s,u;^-) : D^Of, - F{s,u^-,e'J = on Ol, (4.12) 

with terminal condition 0^(7r„;t,x) = 6f^^i{'iTn,{t,x);t,0) on dO^^, where uj'^" : = 
B^Mn'^'^ is deterministic, and 6n{-Kn;T,x) = S.{uj^") when t„ = T. Further, it fol- 
lows from Proposition 7.2 of [6] that for every 6 > 0, there is 6n^ E C^''^{Of^) which 
is a classical supersolution of (4.12) such that On {iTn',t,x) > 0,^_^-^(7r.„, (t, x); t, 0) on 
dOl and \en'^ -0'n\ < S on Ol^. Let 6n = e/2", Hf := H°'°'", and 5^ be the linear 
interpolation of (h|, )«>o- Define 

oo 

i;'{t,uj) := ^(5„ + 0l''5"((Hf,SH|)o<i<n;t,i?t-i?H|,))l[H?,,H^^,), 

n=0 

and denote 

i>-l 

One can check straightforwardly that satisfies the conditions of Definition 4.8 (c), 
■ip'{T,io) > ^{B^), and 

-dti^' -ti{s,B')-d^^P' -a{s,B')■.^^^^P' -F(^s,B',i;'{s,B')^ > 0. (4.13) 
Then Y := ip^, Z := d^ij)'' satisfy the BSDE 

Ys = ^H.+i + F(^r,B'X)dr- Zr(^dBr-li{r,B'yr), ¥x-a.s. 
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on every interval [HjjHj+i) such that supo<t<T \yt — Yt\ < which imphes that (4.10) 
holds true for ip^ and hence ijj'' S Cqq{A^). Notice that \\B^ — B\\t < £, then for some 
constant C, 

\aB')-aB)\ < Ce, \F{s,B',y)-F{s,B,y)\ < Ce. 

Set 

V^:= V^ + 2Ce[l + r-t], 

one can verify that G Cqq{IsP) is a viscosity supersolution of (4.5), and it follows by 
the partial comparison principle in Lemma 4.9 that n^(0, 0) < "0(0,0). Similarly, we 
can construct a viscosity subsolution V S Cq'o(A'^) such that u^(0,0) > '(/'(0,0) and 
IV' ~ V'l ^ 4C[2 + T]e. By sending e — t- 0, we conclude the proof. □ 

4.2 Proof of Theorem 2.12 

Finally, we can complete the proof of our main result which gives a representation of 
BSDE by branching process. 

Proof of Theorem 2.12. By Theorems 4.6 and 4.7, we only need to show that u is a 
viscosity solution of (4.5) with terminal condition ij: and generator F„q defined in (2.14) 
following Definition 4.4. We shall only show the subsolution part. Moreover, we recall 
that in the branching process, the process ^'^X^^^ associated with the first particle is 
extended after its default time Ti by ^'^X^^^ := <^*'''(s, i3*'(^)) for all s G [t,T], where 
is defined by Blf] := W} for all s G [t, T]. 
Suppose that v is not a viscosity subsolution of (4.5), then by Definition 4.4 and 
Remark 4.5, there is (^0,^0) G A*' and G Av{tQ,bjQ) such that v{tQ,ujQ) = (^(to,0) 
and 

- £.'f{tQ,ujo) = - Cip{to,ujo) - l3g(p{to,uJo) = c > 0, (4.14) 
where C is defined by (4.4) and 

no 

Qip{t,oj) := y^^ak{t,oj)ip^{t,uj) - ip{t,oj). 

k=0 

Without loss of generality, we suppose that to = 0. Then '^'•^oxW = o.Ox(i). Further, 
it follows by the continuity of functions 99 and v in Proposition 2.11 that for every e > 0, 
there is n £ such that for every t G [0, f] (with r := h( °'°X^^^)), 

\v{t, O'OX^')) -i;(0,0)| + \gip{t, 0'°X(')) -g^;(t, 0'°X.('))e-'^*| < e. 

Clearly, r is a F*^— stopping time (see (2.11)) in probability space (il,J^, P). Denote 
H/, = H A /i the F°-stopping time on (f]°, Pq), := f A /i and Xt := °'°xf ^ the 
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f''— stopping time and process on (il, J^, P), it follows from equation (3.4) of Lemma 
3.4, together with (4.14), that 

= E^[^(fh,X.)-v{fh,X.) 
= ¥f [ipifh,X.) - ifiO, 0) + v{0, 0) - v{fh,X, 

ip{fh,X.) - ip{0, 0) + gv{Ti,X.)lr,>T, 



{v{n,X.)-v{fH,X.))U^>T, 



< 



E 



(3gip{t,X.) - gv{t,X.)l3e-^^)dt 



+ E^[{v{Ti,X.) - v{fh,X.))lf^>T, 
< {-c + el3 + 2e)E^[fh] < 

for e small enough, which is in contradiction with Definit 
solution. 



ion 4.4 of the viscosity sub- 

□ 



5 Numerical examples 



In this section, we provide a numerical illustration of our representation result, and 
the corresponding numerical implications. Consider the two following PDEs: 

dtvi+xdavi + ^a'^x'^dl^vi+(3{vl-vi) = 0, vi{T,x,a) = i;{x,a) :PDE1 (5.1) 

dtV2 + xdaV2 + ^a^x^dl^V2 + l3{-vl - V2) = 0, V2{T, X, a) = ^{x, a) : PDE2 (5.2) 

with 'ip{x, a) = (a — 1)^. These PDEs correspond to the forward backward SDEs 

dXt = aXtdBt , Xo = 1 (5.3) 
dAt = Xtdt ,^0 = (5.4) 
dYt = -/3{F{Yt)-Yt)dt + ZtdBt , YT = ij{XT,AT) (5.5) 

with the non-linearities Fi{y) = and -^2(2/) = —y^- In our numerical experiments, 
we have taken a diffusion coefficient a = 0.2 and a Poisson intensity (3 = 0.1, and the 
maturity T = 2 or T = 5 years. For T = 2 years (resp. 5 years), the probability of 
default is around 0.18 (resp. 0.39). 

In comparison with the KPP type PDE with F\{y) = y^, the replacement of the 
non-linearity by — has added the term (_)^t-i ^j^g multiplicative functional 
(see Equation (2.12)), without changing the complexity of the branching diffusion 
algorithm. More precisely, we have: 

Nt 

vi{Q,Xo,Aq) = Eo,,.[nV(^^,^^)], (5.6) 

i=l 

Nt 

V2{^,Xo,Ao) = Eo,.[n(-l)^^-V(^^,^^)]- (5.7) 

i=l 
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N 


Fair(PDEl) 


Stdev(PDEl) 


Fair(PDE2) 


Stdev(PDE2) 


CPU (seconds) 


12 


5.69 


0.16 


5.36 


0.16 


0.1 


14 


5.61 


0.08 


5.23 


0.08 


0.6 


16 


5.50 


0.04 


5.15 


0.04 


1.5 


18 


5.52 


0.02 


5.16 


0.02 


5.9 


20 


5.53 


0.01 


5.16 


0.01 


23.6 


22 


5.54 


0.00 


5.17 


0.01 


94.1 



Table 1: MC price quoted in percent as a function of the number of MC paths 2 . PDE 
pricer(PDEl) = 5.54. PDE pricer(PDE2) = 5.17 (CPU PDE: 10 seconds). Maturity= 
2 years. Non-hnearities for PDEl (resp. PDE2) Fi{u) = v? (resp. F2{u) = —u^). For 
completeness, the price with /5 = (which can be obtained using a classical Monte-Carlo 
pricer) is 6.52. 

Our branching diffusion algorithm has been checked against a two-dimensional PDE 
solver with an ADI scheme (see Tables 1,2). The degenerate PDEs have been converted 
into elliptic PDEs by introducing the process At = f^Xgds + {T — t)Xt, satisfying 
dAf = (T — t)dXt. The computational experiments was done using a PC with 2.99 
Ghz Intel Core 2 Duo CPU. 

Note that our algorithm converges to the exact PDE result as expected and the error 
is properly indicated by the Monte-Carlo standard deviation estimator (see column 
Stdev). In order to illustrate the impact of the non-linearity F on the price f , we have 
indicated the price corresponding to /? = 0. 

We would like to highlight that replacing the generator |x^o"^9^ -|- xda by a multi- 
dimensional operator C can be easily handled in our framework by simulating the 
branching particles with a diffusion process associated to C. This is out-of-reach with 
finite-difference scheme methods and not such an easy step for the BSDE approach 
which require computing conditional expectations. 
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0.03 
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22 
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